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ABSTRACT 

Recent studies have presented evidence that the Milky Way global potential may be non- 
spherical. In this case, the assembling process of the Galaxy may have left long lasting stellar 
halo kinematic fossils due to the shape of the dark matter halo, potentially originated by orbital 
resonances. We further investigate such possibility, considering now potential models further 
away from ACDM halos, like scalar field dark matter halos, MOND, and including several other 
factors that may mimic the emergence and permanence of kinematic groups, such as, a spherical 
and triaxial halo with an embedded disk potential. We find that regardless of the density profile 
(DM nature), kinematic groups only appear in the presence of a triaxial halo potential. For the 
case of a MOND like gravity theory no kinematic structure is present. We conclude that the 
detection of these kinematic stellar groups could confirm the predicted triaxiality of dark halos 
in cosmological galaxy formation scenarios. 

Subject headings: Galaxy: halo — Galaxy: kinematics and dynamics — Galaxy: structure — 


1. Introduction 


The ACDM scenario is considered the stan¬ 
dard one in cosmology because it incorporates self- 
consistently many large scales properties of the 
universe. An intrinsic feature of this scenario is 
the hypothesis of the existence of cold dark mat¬ 
ter particles required to explain the cosmic struc¬ 
ture formation. The ACDM model is also able 
to make predictions on galactic scales, relating 
the process of galaxy formation with the internal 
properties of galaxies, triggering a vivid debate in 


the astronomical community, (e.g. [White fc Rees 


Bosch|1998 Ghigna et al.|1998 

Mypin et al.|1999t 

Moore et al. 1999 

Avila-Reese 

20061 Pizagno et 

al. 2007| jValenzue 

a et al.||2007 

Covernato et al. 


2012 ). 


Cosmological simulations had shown that the 


formation of galactic size dark matter halos occur 
through the hierarchical assembly of CDM struc¬ 
tures ( Kauffmann et al.||1993 Klypin et al.||1999 


Moore et al.|1999 ), where small systems form first 


and then merge together to form more massive 
ones. This process leads to halos with a strongly 


triaxial shape (Allgood et al. 2006 

ir] [MT] ). 


Vera-Giro et 


Therefore, the possible detection of the dark 
matter halos triaxiality would be of great impor¬ 
tance for these galactic evolution and formation 
theories, as it would confirm one of their most in¬ 
trinsic predictions. 


In paper I (Rojas-Niho et al. 2012), we intro¬ 


duced a new strategy in order to detect triaxiality 
in the Milky Way halo, in addition to the ones 


commonly discussed in literature (Jing & Suto 
|2002[ [Law et ar]|2009| jPeharrnbia et al.||2009| 
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1 IVera-Ciro 20141 Deg & Widrow 20131 Lux 

et al 

.. 12012). In a more recent paper Loebman et 

al. ( 

2014) solve the Jeans equation to constraint 


the shape and distribution of the dark matter halo 
within the Milky Way and find out that the solu¬ 
tion is consistent with an oblate halo. The former 
method is very promising, although it has an im¬ 
portant dependence on assumptions currently val¬ 
idated mostly by cosmological MW formation sim¬ 
ulations which are still uncertain. Our work is in¬ 
dependent and complementary to all of them. Our 
strategy is based on the stellar kinematics in the 
Galaxy halo. The idea was to compare the orbital 
structure generated by a spherical halo with the 
orbital structure generated with a triaxial halo. 
The kinematic differences between both cases and 
their comparison with observational data could 
help to determine whether the Galactic halo is re¬ 
ally triaxial. An important difference between a 
triaxial potential and a spherical one is that in the 
first case there are abundant resonant orbits, while 
in the second rosette orbits are dominant and have 
no angular preference. The presence of resonant 
orbits favour the formation of moving groups of 
stars: the kinematic groups. 

In paper I we showed that if the dark matter 
halo of the Milky Way possesses a triaxial shape, 
the assembling history of the stellar halo is able to 
populate the quasi-resonant orbits triggering kine¬ 
matic stellar groups. 

The analysis of the resulting kinematic and or¬ 
bital structure may give evidence for the halo tri- 
axiality. 


For the study presented in paper I we only used 
the gravitational potential of a NFW dark mat¬ 
ter halo with a density profile that results from 
CDM-only simulations ( Navarro et al.||1996 ). But 
regardless of the great success of the ACDM sce¬ 
nario explaining the large scale structure of the 
universe, the lack of clear positive results in at¬ 
tempts of direct dark matter detection have mo¬ 
tivated the exploration of alternative models. In 
this way, with the purpose of testing other theo¬ 
ries, alternative to cold dark matter and their ef¬ 
fect on the dynamical structure of the Galaxy, we 
introduce two different models for the potential, 
one based on modified gravity theories (MOND 


specifically), and a different halo model motivated 
on Scalar Field Dark Matter (SFDM). 


Regarding to MOND, the work of IMilgronT"! 
(1983) is one of the most popular proposals in 
astrophysics. A common feature in these fam¬ 
ily of theories is that the only gravity source are 
baryons. In the Milky Way case, that implies that 
mostly, the flat disk potential is responsible of the 
whole galactic dynamics, particularly in the stellar 
halo. 


On the other hand, regarding the SFDM, this 
is an alternative scenario that has received much 
attention recently, the main hypothesis is that 
the nature of the dark matter is determined by a 
fundamental scalar field that condensates form¬ 
ing Bose-Einstein Condensate (BEG) “drops” 


(Guzman & Matos 2000 Magana et al.| 2012), 


these condensates represent galaxies dark matter 
haloes. Robles & Matos (2013) consider that dark 
matter is a self-interacting real scalar field em¬ 
bedded in a thermal bath at temperature T with 
an initial Z 2 symmetric potential. Due to the ex¬ 
pansion of the Universe, the temperature drops 
in such a way, that the Z 2 symmetry is sponta¬ 
neously broken, and the field rolls down to a new 
minimum. 


To study the galactic scale of the model one 
can solve the Newtonian limit of the equation that 
describes a scalar field perturbation, i.e., where 
the field is near the minimum of the potential 
that describes its interaction and where the grav¬ 
itational potential is locally homogeneous. Eor 
galaxies, this Newtonian approximation provides 
a good description giving an exact analytic so¬ 
lution within the Newtonian limit. The follow¬ 
ing density profile represents halos in condensed 
state or halos in a combination of exited states, 
p = pQ {kr)/(krY ^ where po and k are fit¬ 
ting parameters. Efforts to include the baryonic 
influence on SEDM halos are necessary to have 
accurate comparisons with observations however 
they are still in very early stages of understanding 
( Gonzalez-Morales et al.|[2M3 ). 

An important feature of SEDM density profiles 
is the presence of wiggles, characteristic oscilla¬ 
tions of scalar field configurations in excited states. 
This could result in some differences in the gravi¬ 
tational potential for a scalar field halo compared 
with the ACDM profile, possibly imprinting a sig¬ 
nature in the stellar kinematics of the galactic stel- 
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lar halo. 

The aim of this paper is to generalize the results 
obtained in paper I, where we only employed tri- 
axial and spherical NFW halos. Now we carry out 
numerical simulations with other galactic compo¬ 
nents as a Miyamoto-Nagai disk and a different 
dark matter density profile, and other gravity the¬ 
ory (MOND disk). Finally, for this work we also 
generalize the initial conditions used in paper L 

This paper is organized as follows. In Section 
the 3-D galactic potentials used to compute orbits 
are briefly described. In Section we introduce 
our numerical simulations, techniques and a strat¬ 
egy aimed to efficiently explore the stellar phase 
space accessible to a hypothetical observer, we also 
present the results of our numerical simulations. 
In Section we generalize our initial conditions 
scheme. Finally, in Section we present a discus¬ 
sion of our results and our conclusions. 


tern. However we emphasize that the study 
does not need a detailed model of the Milky 
Way Galaxy. 


2.1. Halo Potentials 


We have implemented two intrinsically different 
potential halos, one as in Paper I, produced by a 
NFW density profile ( Peharrubia et al.]|2QQ9 ) 


^{x^y^z) = 27rGabcporl 

X r (1) 

Jo + s(t) y(a2 +t)(&2 + t)(c 2 + r)’ 

where po is the characteristic halo density, the di¬ 
mensionless quantities a, b and c are the three 
main axes, and Vg is the radial scale. The tri- 
axiality effect is accounted by using elliptical co¬ 
ordinates, where. 


2. The Potential Models 


In Paper I we assumed a steady NFW (Navarro 
et al.||1996 ) triaxial dark matter halo as a proof of 


concept that a triaxial halo develops and preserves 
abundant structure in the stellar phase space be¬ 
cause of the resonant orbital structure. It is worth 
mentioning that this is applicable also to other 
non-spherical halo profiles (oblate, prolate). With 
that study we were able to produce clear features 
in the velocity space (i.e. halo moving groups), 
but we did not include any other galactic compo¬ 
nents or halo representations other than a galactic 
triaxial NFW profile halo. 


In this paper, with the purpose of testing the 
generality of the results in Paper I we have ex¬ 
tended our studies including now the effect of a 
disk and a very different triaxial halo potential. 
We have also explored a disk that responds to a 
modified law of gravity (MOND). 

Our results are applicable regardless the 
specific values of the gravitational potential 
parameters, they only depend on the triax- 
iality, therefore the kinematic signature is 
expected either in the Milky Way or other 
galaxies that may have non-spherical halos. 
As the best known case is the Milky Way, 
and the one that nearest in the future will 
have stellar kinematics observations avail¬ 
able, we will consider it as our testbed sys- 


s(t) 


+ T + T 



( 2 ) 


For the NFW halo the density profile depends on 
the free parameters r^, po and the axis ratios a, b 
and c. As in Paper I, the rotation curve is used 
as one of the primary observational constraints 
for the parameters, so we kept Vg = 8.5 kpc, and 
Po = 0.056M© obtained assuming a maxi¬ 

mum rotation velocity of 220kms“^. 


For the second halo we used the density pro¬ 


file of a scalar field configuration (Robles & Matos 


2013) and included the triaxiality using the for¬ 
mulae of Chandrasekhar (1969), that leads to the 
next gravitational potential 


^{x,y,z) = 

irGabcpo ln{s{T)) — Gi{2ks{r)) 

Jo + T^/b‘^ + + t 

where a, b, and c are the three main axes and Gi 
is the cosine integral function. 

For the scalar field halo we chose the parame¬ 
ters k = 0.123 kpc“^ and po = 0.0196 Mq pc~^ 
to guarantee a maximum rotation velocity of 220 
kms“^. It is worth mention here that, with 
a correct selection of the parameters, the 
SFDM model has proved to be in a good 
agreement with observational data at galac¬ 
tic scales, particularly with rotation curves. 
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And for several cases, a SFDM halo fits even 
better data than a NFW halo (see e.g. figs. 
5,6,7 and 8 of | Martmez-Medina Matos | 
(( 20 ^). 

2.2. Disk Potential 

With the halo potential, we have included a 
potential that simulates the Galactic disk. We 
considered two cases, a Miyamoto-Nagai and a 
Kuzmin model for a modified gravity case. With 
these models we performed numerical simulations 
in order to explore the influence of the disk on the 
orbital structure of the Galaxy. 


2.2.1. Miyamoto-Nagai potential 


A commonly used model for the Galactic disk 
is the one proposed by Miyamoto & Nagai (1975), 


^R,z) = - 


GM 


+ (a + + 6^)^ 


( 4 ) 


This potential has three free parameters, a and b 
that represent the radial and vertical length-scale 
respectively, and M the total mass of the disk. For 
the case of the Milky Way we adopted a = 5.3178 
kpc, b = 0.25 kpc, and M = 8.56 x 10^^ M0 
( Allen fc Santillan~||1991 ), values that are in good 
agreement with observations. 

2.2.2. MOND Kuzmin disk potential 

With the purpose of searching for a difference 
in the orbital structure (moving groups) induced 
by using the typical models for the Galactic po¬ 
tential, that assume the existence of dark mat¬ 
ter massive triaxial halos and modified newtonian 
dynamical models for the gravity, we have con¬ 
structed a MOND (Modified Newtonian Dynam¬ 
ics) galactic disk as proposed by Milgrom (1983), 
as an alternative to the dark matter halo. 

Unlike newtonian gravity, the MOND version of 
Poisson equation is non-linear and extremely dif¬ 
ficult to solve. However, for some mass distribu¬ 
tions, it is possible to find the gravitational poten¬ 
tial and a Kuzmin disk is one of these cases. The 
advantage of employing a Kuzmin disk is that the 
MOND gravitational field can be obtained from 
the newtonian field ( Read fc Moore|[2QQ5 ), 


g = 9 n 




l + {l+Aal/glY/^ 
2 


( 5 ) 


where ao is the MOND constant, qn the newto¬ 
nian gravitational field generated by the Kuzmin 
disk, and g is the corresponding MOND gravita¬ 
tional field. 

Gonsequently, since the Miyamoto-Nagai po¬ 
tential reduces to a Kuzmin potential when b=0, 
that corresponds to the case of a completely flat 
disk. The potential produced in this case is 


^(R^z) =- , =. 

yi?2 + (a+|z|)2 

Thus is calculated first from the Kuzmin po¬ 
tential, and then we obtain g with equation (§. 
Once we have we carry out the integration. 

In the next section, we explore the effect of the 
halo and the disk potential in the resulting orbital 
structure. 


Numerical Test Particle Simulations 


Our physical system consists of test particles 
moving under the influence of a fixed gravitational 
potential which is the superposition of different 
galactic components as described in the sections 
above. The temporal evolution of the system is 


solved using the Bulirsh-Stoer method (Press et 


al. 1992), that works best when following motion 


in smooth gravitational fields. Using this tool we 
follow, for each particle, the evolution of x(t), y{t) 
and z{t) by integrating the equations of motion nu¬ 
merically. As shown in paper I, the orbital analysis 
is performed in the velocity space focusing in the 
kinematic projections Vx - by means of an or¬ 
bital classification. For this purpose we apply the 


spectral method by Garpintero & Aguilar (1998) 


that automatically classifies an orbit, distinguish¬ 
ing it between a regular and an irregular one. It 
can also identify loop, box, and other resonant or¬ 
bits, as well as high order resonances. 

In our previous work we already saw that for a 
triaxial halo, resonances trigger kinematic stellar 
structure in the case of a single accretion event and 
also when the potential is populated in a stochas¬ 
tic way trying to mimic the stellar halo assembly 
history. 
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In order to generalize the results seen in paper 
I, we follow two approaches to generate the initial 
set ups. The first one is the one already used in the 
previous paper, where we choose in advance, the 
observation point. The particles are placed ran¬ 
domly around this point with velocities ranging 
between zero and the local escape velocity. This 
set up will mimic the contribution of many accre¬ 
tion events with different orientations and energies 
but keeping only bounded orbits. 

This method has proved to be efficient to ex¬ 
plore the available phase space, as it will be con¬ 
firmed in the next sections when using more gen¬ 
eral initial conditions. This is because a star that 
falls in a resonance will return to the neighbour¬ 
hood of the observation point, therefore, we are 
focusing on orbits that the artificial observer will 
detect. The process produce persisting kinematic 
groups at the velocities around resonances because 
irregular and open orbits will almost never return. 
In contrast a non-resonant region will present 
nearly evenly distributed types of orbits as a result 
of phase mixing, and also non-kinematic groups 
after some mixing time scales. This method in¬ 
creases our particle statistics without the use of 
high computational times. 


3.1. Triaxial halo with disk 


The observation neighbourhood is defined as a 
sphere of 1 kpc of radius and centered at the given 
observation point, in this case {x = 15 kpc, ^ = 0, 
2 ; = 0). We set 2 x 10^ stars randomly distributed 
inside this sphere with velocities that range be¬ 
tween zero and the escape velocity. Because we are 
focusing in the stellar halo, stars generated with 
Vz ^ 0 {vzivx < tan 10°) are removed from the ini¬ 
tial conditions. These stars will move close to the 
X-Y plane, remaining in this plane all along the 
simulation and will always have Vz close to zero. If 
we have populated the whole halo this orbits con¬ 
centration would be compensated with neighbour 
orbits thrown out from different galaxy positions. 
This orbits would appear in the velocity plane as 
an horizontal band but they do not form a kine¬ 
matic group, as we verified it with the spectral or¬ 
bit classifier ( Carpintero fc Aguilar||1998 ). These 
particles are an artifact, already explained in a 
previous letter ( [Rojas- Nino et al. 2012) and should 
not be confused with the ones induced by the pres¬ 
ence of the disk (see Fig. [^. The system evolves 


for some 12 Gigayears and at the end of this in¬ 
tegration time, we studied the stellar kinematic 
distribution. This integration time is long enough 
to be comparable to the Milky Way halo evolu¬ 
tion ( Kalirai||2Q12 ). For our purpose we consider 
the stars that are found after the simulation inside 
the “artificial observer neighbourhood”. In Figure 
[ 1 ] we show the kinematic distribution of stars at 
the end of the simulation for two different cases as 
seen by an observer at 15 kpc from the halo cen¬ 
ter. The upper panel of Figure corresponds to 
a spherical NFW halo (a = 1.0,6 = 1.0, c = 1.0), 
the bottom panel corresponds to a triaxial NFW 
halo, where in order to generalize the results of pa¬ 
per I, we adopt the same values for the axis ratios 
(a = 1.47, b = 1.22, c = 0.98), consistent with cos¬ 
mological simulations (Jing & Suto (2002) 

Giro et al.| ( |MT] )). 


Vera- 


Both cases contain a Miyamoto-Nagai disk as 
described in section 2.2. 1| The Miyamoto-Nagai 
disk is axisymmetric, which is not necessarily true 
in the presence of a triaxial halo, since the poten¬ 
tial of the latter could deform the disk. However, 
this does not affect the results because an ellipti¬ 
cal disk aligned with the major axis of the halo 
only reinforce the kinematic signature of triaxial- 
ity. Note that the only structure that appears in 
the first case is a band close to Vz = 0, corre¬ 
sponding to the stars moving in the X-Y plane. 
As explained in Paper I, although these orbits 
would appear in the velocity plane as a horizontal 
band, they do not form a kinematic group. In the 
second one we can observe symmetrical bands at 
Vz ^ ±80 kms“^. This figure shows a clear differ¬ 
ence between both cases. The symmetrical bands, 
associated to resonant orbits, would be evidence 
for the presence of a triaxial halo. 


3.1.1. Changing the axes. 

The previous axes choice for the triaxial dark 
matter halo allows us to compare with the results 
of paper I. 

Now we change the axes ratios (a = 1.47, 

b = 0.97, c = 0.48), and orientate the disk and 
the dark matter halo according to the halo model 


that fits better the Sagittarius stream (Law et al 


2009). In this way the minor and major axes lie in 


the galactic plane with the intermediate axis per¬ 
pendicular to the disk and with the new val¬ 
ues af the axes we obtain, hp/ap = 0.83 and 
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Cp/ttp = 0.67, corresponding to the axis ratios of 
the potential. 
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Fig. 1.— Stellar kinematic distribution in the 
velocity space, as seen by an observer at 15kpc 
from the halo center, for a spherical NFW halo + 
Miyamoto-Nagai disk (top panel), and for a tri- 
axial NFW halo + Miyamoto-Nagai disk (bottom 
panel). 


Figure shows the result of this experiment, 
with an orientation of the axes different from the 
treated above but motivated by observational con¬ 
straints. The track of the halo triaxiality is even 
more visible for this particular orientation as now 
four horizontal bands develop in the kinematic 
stellar structure. Applying the stellar orbit classi¬ 
fier to one of this bands, that with « 50 km s“^, 
we can establish a resonant origin for these struc¬ 
tures. 

Comparing Figures [^(bottom panel) and[^ the 
differences in the kinematics induced by the two 
configurations is noticeable. Both, however con¬ 
tain resonant orbits in spite the different orienta¬ 
tions of the halo relative to the disk. 

3.2. Kuzmin disk with MOND 

In this section we describe the numerical sim¬ 
ulations carried out considering a Kuzmin disk 


within MOND theory (Milgrom 1983) as an al¬ 


ternative to dark matter. 

In this theory, dark matter does not exist at all. 
Instead it is proposed that Newtonian dynamics 
should be modified to explain observations. In the 
MOND version, gravity law has a more general 
form, 


g GM 
ao 


( 7 ) 


where ao is a universal constant that determines 
the transition between the regime of strong and 
weak field and its value is approximately 1.2 x 
IQ-io ^ g- 2 ^ function g{x) is not determined 
in this theory, but it must satisfy the condition 
jii{x) ~ 1 when x ^ 1 and /i(x) ~ x when x <C 1. 

In the outer part of the galaxies, the weak field 
regime is valid and g{x) ^ x. This leads to a 
rotation velocity v = (GMuq)^^^, therefore, the 
rotation velocity in the outer part of galaxies is in¬ 
dependent of r. Thus MOND can explain the flat¬ 
tening of the rotation curve at great distances from 
the center of the galaxies without invoking the ex¬ 
istence of dark matter. This is its main achieve¬ 
ment and one of the the reasons it was born. 
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However MOND also has several problems, 
especially to correctly describe the dynamics of 
galaxy clusters, for example ( Sander s]| 1999 ; Clowe 
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Fig. 2.— Vx-Vy projection of the velocity space as 
seen by a hypothetical observer placed at 15kpc 
from the center in the case of a NFW triaxial halo 
+ Miyamoto-Nagai disk with the axes orienta- 

(2009). The histogram 


tion taken from Law et al. 


shows the result of applying a spectral orbital clas¬ 
sifier were we see that the regions more densely 
populated correspond to resonances. This sup¬ 
ports the interpretation that the kinematic struc¬ 
ture has a resonant origin. The upper inset shows 
the orbit of a single particle inside one of these 
resonant regions. 


Nevertheless, despite that dark mat¬ 
ter is the most accepted theory, MOND (or mod¬ 
ifications of it), has not been totally ruled out. 

For this purpose we construct initial conditions 
in the same way as in the cases of the triaxial halo 
and a disk: 2 x 10^ stars within a sphere of 1 kpc 
radius but with the observer neighbourhood cen¬ 
tered at {x = 15 kpc, ^ = 0, 2 : = 0), with random 
velocities between zero and the escape velocity. 

We analyze the resulting velocity space after an 
integration time of 12 gigayears. Figure shows 
the kinematic structure for this case. An area 
of higher density in velocity space is formed, but 
without the symmetric bands that characterize the 
case of a triaxial halo. 

As shown also in Figure a classification of or¬ 
bits give us more information about the kinematic 
structure, with no presence of stellar groups and 
dominated by loop orbits. 

Although the geometry of the kinematic struc¬ 
ture is preserved regardless of the detailed model 
of the disk, nevertheless, we decided to repeat the 
calculations with a Kuzmin disk within a triaxial 
halo with newtonian gravity and found that the 
structure is very similar to the one shown in Fig¬ 
ure where we used a Miyamoto-Nagai disk, we 
show this result in Figure]^ 

4. Generalized Initial Conditions 

As described above the initial conditions set¬ 
up, used in the previous sections and in paper I, 
provide an observer, with a good kinematic accu¬ 
racy, only inside a limited vicinity, and avoiding to 
populate the whole halo with particles, ensures a 
relatively low computing cost. However, this ini¬ 
tial distribution of particles is not isotropic around 
the halo center and the question remains about 
whether this choice for the initial conditions has a 
non-negligible effect on the final kinematic struc¬ 
ture for the triaxial dark matter halo employed in 
paper I. With this in mind we introduced a more 
general and complete initial conditions set-up, dif¬ 
ferent to the described above. 

In order to do this, we consider initial condi¬ 
tions with spherical symmetry around the halo 
center populating the entire simulated domain. 
This is done by increasing the number of particles 
which are distributed inside a sphere of 25 kpc of 
radius centered at (x = 0, ^ = 0, 2 : = 0). The ve- 
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Fig. 3.— Vx—Vz projection of the stellar kinematic 
distribution for a Kuzmin disk within MOND the¬ 
ory. The histogram shows the result of the orbit 
classification. The upper inset shows a loop or¬ 
bit in its projections x — y and x — z for a single 
particle inside the classified region. 


locities were selected randomly between zero and 
the local escape velocity. 

Although computationally expensive, sampling 
the entire halo and increasing the number of par¬ 
ticles give us a more conclusive mapping of the 
gravitational potential to establish its relationship 
with the distribution of stars in velocity space. 
The choice for the initial conditions ensures that 
the kinematic stellar structure is not due to the 
sampling but due to the physical system. 

With these new initial conditions, we now re¬ 
peat the experiment of paper I with just a triaxial 
NFW halo. Finally, a different density profile for 
the dark matter halo motivated within the scalar 
field dark matter model is also employed. The 
generation of stellar groups within these two, very 
different in nature, triaxial halo potentials allow 
us to establish its origin in the triaxiality of the 
halo rather than in the dark matter nature. 

4.1. Triaxial NFW halo 

In this experiment we use the potential of equa¬ 
tion Q with parameters (a = 1.47,6 = 1.22, c = 
0.98) for the axis of the triaxial halo. 

We used the generalized initial conditions de¬ 
scribed above, 16 x 10^ particles populating a 
sphere of 25kpc of radius centered at {x = O^y = 
O^z = 0) with velocities selected randomly be¬ 
tween zero and the local escape velocity. The sys¬ 
tem evolves 12 Gyr. 

Figure shows the Vx-Vz projection of the ve¬ 
locity space at the end of the integration time 
as seen by a hypothetical observer placed at 
15kpc from the halo center. Overdensities are 
distinguishable in the kinematic structure, mostly 
around Vz ^ 45kms“^ and Vz ^ —45kms“^. 
By doing a classification of orbits one can notice 
that these overdensities have a resonant origin and 
its independent of the initial conditions. 

4.2. Triaxial scalar field dark matter halo 

For this triaxial halo (equation we used 
the generalized initial conditions described above, 
12 X 10^ particles populating a sphere of 25kpc 
of radius centered at (x = 0, ^ = 0, 2 ; = 0), with 
velocities selected randomly between zero and the 
local escape velocity. 

We selected (a = 1.47,6 = 1.22, c = 0.98) for 
the triaxial halo axes, and evolved the system for 
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12 Gyr. Figure shows the Vx-Vz projection of 
the velocity space at the end of the integration 
time, as seen by a hypothetical observer placed at 
8.5 kpc from the halo center. 


In figure 1^ we found horizontal structure resem¬ 
bling that seen in our previous experiments. We 
also performed a simulation with the spherical case 
{a = 1,5 = l,c = 1) for comparison, and found 
that the velocity space is featureless and homoge¬ 
neous. This suggests that the kinematic structure 
seen in the velocity space of the triaxial halo corre¬ 
sponds, again, to stars librating in the vicinity of 
resonances. We verified this by using the spectral 
orbital method of Carpintero & Aguilar (1998), 
finding that the most remarkable regions of the 
velocity space are dominated by resonant orbits. 


In figure we show the Vx-Vy projection for 
another measure of the data but now moving the 
hypothetical observer to a distance of 20 kpc from 
the halo center. At this distance we found again 
a horizontal pattern as in the figures before which 
means that the effect of triaxiality acts across the 
entire stellar halo. Again by applying the spectral 
orbital method we can see that dark matter halo 
resonances are able to densely populate the stellar 
velocity space with this characteristic horizontal 
pattern. 


5. Discussion and Conclusions 


In this work we extend the orbital study of a 
stellar system under the influence of the gravita¬ 
tional potential generated by a steady triaxial dark 
matter halo, compared to spherical ones as in Pa¬ 


per I (Rojas-Niho et al. 2012). We now have added 


a Milky-Way-like disk to our study in order to in¬ 
vestigate for a way to separate the effects from a 
non-spherical halo and a disk. We have also tested 
a different type of dark matter halo profile based 
on theories of scalar field dark matter. Finally, we 
have included a Kuzmin disk assuming a MOND 
gravity theory in order to test the consequences on 
the orbital behaviour and the possible production 
of moving groups. 

We show that, independently of the nature of 
dark matter, a non-spherical shape of the Milky 
Way dark matter halo, influences strongly the stel¬ 
lar halo kinematic structure, as seen from the long 
lasting resonant features formed in simulations. 
As it was pointed out in Paper I, there is an enor¬ 


mous difference with the spherical case, where no 
kinematic structure is triggered. We have explored 
also with a different halo profile motivated by the 
scalar field dark matter model with an important 
feature, the presence of density wiggles, character¬ 
istic oscillations of scalar field configurations. We 
conclude that this scalar field density wiggles are 
unable to trigger any moving groups in the spher¬ 
ical case. However in the triaxial case, this profile 
produces structures with the same resonant origin 
as in the typical NFW triaxial halos. 

We also confirm that the structures seen in the 
velocity space for triaxial halos, do not depend 
on the initial conditions, the effect is still present 
in spite of the two different initial distributions 
strategies of particles used in this paper. The fact 
that the kinematic structure is independent from 
the initial conditions makes stronger the resonant 
nature versus an incomplete relaxation of initial 
conditions. Our ICs strategy is motivated by the 
fact that there are not selfconsistent phase space 
distribution for triaxial halos, it can be argued 
that incomplete relaxation is producing the kine¬ 
matic structure reported in this paper. However 
the fact that such structure is produced only in 
the non-spherical case, and lasts for longer than 
the Universe age, and is not altered by the stellar 
disk presence makes our conclusions robust. 

In this work a disk potential was introduced 
in order to understand whether this component 
is able to erase the orbital structure produced by 
a triaxial halo, or if the disk itself could produce 
moving groups in the halo that can be erroneously 
taken as the ones produced by a triaxial halo. 
With this study we learnt that the disk influence 
is only important close to it, and that it is unable 
to erase the orbital structure induced by the halo 
in all cases we studied. In the same manner, we 
studied the effect of the disk on the production 
of kinematic structures in the halo, and we found 
that for an axisymmetric disk, there are no dom¬ 
inant resonant structures (moving groups). Al¬ 
though it is known that non-axisymmetric struc¬ 
tures (such as spiral arms and bars) in the disk, 
may produce important kinematic structure into 
galactic disks, their influence diminishes fast out 
of the disk plane. 

Finally, we have also tested the results compar¬ 
ing with a modified gravity model (MOND), with 
a Kuzmin disk. We find in this case no kinematic 
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structure at all. 

In the close future we will be likely able to dis¬ 
tinguish, for instance, satellite remnants groups 
originated from a single event from these associ¬ 
ated with the dark matter halo shape with the 
help of the stellar population data. The great sur¬ 
veys such as Gaia, SDSS, etc., will allow us the 
detection of these type of kinematic groups, which 
would be a clear indication of the triaxiality (non¬ 
sphericity) of the Milky Way dark matter halo, 
and will be of great importance for the galaxy for¬ 
mation theories and the ACDM scenario. 
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Fig. 4.— Vx—Vz projection of the stellar kinematic 
distribution for a triaxial halo with a Kuzmin disk 
within Newtonian gravity. The presence of kine¬ 
matic structures is clear in this plot. 


This 2-column preprint was prepared with the AAS IATeX 
macros v5.2. 
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Fig. 5.— Vx-Vz projection of the velocity space 
as seen by a hypothetical observer in the case of 
a NFW triaxial halo. After evolving for about 
12 Gyr, we found kinematic structure as seen in 
paper I, even with the different initial conditions. 
The histograms show the result of applying a spec¬ 
tral orbital method were we see that the regions 
more densely populated correspond to resonances. 
This supports the interpretation that the kine¬ 
matic structure has a resonant origin. The upper 
inset shows the orbit of a single particle inside one 
of these resonant regions. 


Fig. 6.— Vx-Vz projection of the velocity space as 
seen by a hypothetical observer in the case of the 
scalar field triaxial halo. After evolving for about 
12 Gyr, we found kinematic structure as seen in 
paper I, even with the randomly selected veloci¬ 
ties and homogeneous initial particle distribution 
around the halo center. The histograms showed 
here are the result of applying a spectral orbital 
method, where regions more densely populated, 
correspond to resonances. The upper inset shows 
the orbit of a single particle inside one of these 
resonant regions. 
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Fig. 7.— Vx-Vy projection of the velocity space 
as seen by a hypothetical observer in the case of 
the scalar field triaxial halo. Taken from the same 
simulation as before but this time the observer is 
placed at a distance of 20 kpc from the halo center. 
The histograms showed are the result of applying 
a spectral orbital method were we see that the re¬ 
gions more densely populated correspond to reso¬ 
nances. This supports the interpretation that the 
kinematic structure has a resonant origin. The 
bottom inset shows the orbit of a single particle 
inside one of these resonant regions. 
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